Dynamics of ultracold molecules in confined geometry and electric field 
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We present a time-independent quantum formalism to describe the dynamics of molecules 
with permanent electric dipole moments in a two-dimensional confined geometry such as a one- 
dimensional optical lattice, in the presence of an electric field. Bose/Fermi statistics and selection 
rules play a crucial role in the dynamics. As examples, we compare the dynamics of confined 
fermionic and bosonic polar KRb molecules under different confinements and electric fields. We 
show how chemical reactions can be suppressed, either by a "statistical suppression" which applies 
for fermions at small electric fields and confinements, or by a "potential energy suppression" , which 
applies for both fermions and bosons at high electric fields and confinements. We also explore col- 
lisions that transfer molecules from one state of the confining potential to another. Although these 
collisions can be significant, we show that they do not play a role in the loss of the total number of 
molecules in the gas. 



I. INTRODUCTION 

Experimental evidence for ultracold chemistry of 
quantum-state controlled molecules [l[ and dipolar col- 
lisions in the quantum regime has been obtained re- 
cently for fermionic KRb molecules in the lowest elec- 
tronic, vibrational, rotational quantum state Q and well- 
defined hyperfine states 3 . Bosonic species of KRb have 
also been formed recently H as well as other alkali polar 
molecules such as RbCs Q and Li Cs Ml • The exoergic re- 
action KRb + KRb -> K 2 + Rb 2 0413 prevents long trap 
lifetimes of these molecules, especially in electric fields, 
where the chemical reactivity increases as the sixth power 
of the dipole moment induced by the electric field 0, [H| ■ 
Lifetimes are then typically of the order of 10 ms for ex- 
perimental electric fields. However, polar molecules offer 
long range and anisotropic dipolar interactions in elec- 
tric fields. If the molecules are confined in optical lat- 
tices, they can be stabilize against collisions and chem- 
ical reactions [H|-[I3, if the dipoles are polarized in the 
direction of a tight confinement. If these molecules are 
confined into the ground state of a realistic one dimen- 
sional optical lattice, electric field suppression of chem- 
ical reactions is expected to occur, yielding lifetimes of 
KRb molecules of ~ 1 s and elastic scattering rates 100 
times more efficient than chemical reaction rates [HI, [lj| . 
Both of these are needed to achieve molecular evapo- 
rative cooling and to reach the quantum regime where 
the phase-space density is high. For fermionic molecules, 
creation of degenerate Fermi gases of dipoles will likely 
be possible. In case of bosonic molecules, Bose-Einstein 
condensates can instead be formed. This will reveal ex- 
citing physics with ultracold controlled molecules in the 
quantum regime fl~8l - l2l| . 

We address in this paper two important points regard- 
ing collisions in a lattice. First, suppression of confined 
chemical reactions in electric fields can be obtained by 
using the centrifugal repulsion of fermionic molecules 
in the same internal state (electronic, vibrational, ro- 
tational and spin) and in the same confining state of 
the one dimensional optical lattice. The centrifugal re- 



pulsion comes from the statistics of identical fermions 
in indistinguishable states. This requires only compara- 
tively small dipoles and weak confinements. Suppression 
that relies directly on the confining potential and the 
repulsion due to electric dipoles can also be obtained, 
but requires larger dipoles and stronger confinements. It 
does, however, suppress both bosons and fermions, in in- 
distinguishable states or not, or even for different polar 
molecules. 

Secondly, realistic experimental dynamics of polar 
molecules in confined geometry is more complicated than 
the ideal case used in the recent theoretical works [l5l . [l6| . 
where only molecules in the ground state of the lattice 
were considered. Realistically, molecules can also be 
formed in excited states of the optical lattice, depend- 
ing for example on the temperature, the strength of the 
confinement, and the way the optical lattice is turned 
on 22] . It is therefore important to know how rapidly col- 
lisions can populate higher confining states, which could 
after all, contribute to re-thcrmalization; and how do the 
molecules in these excited states affect the loss rate of 
the total molecules. These questions are important for 
ongo ing experiments of KRb molecules in an optial lat- 
tice |22| |. 

In this article, we extend the formalism developpcd 
in our former work [lj|. We describe in section II the 
dynamics of molecules in an arbitrary initial confining 
state of the lattice, and consider the possibility for the 
molecules to leave such a state for another after a colli- 
sion. In section III, we show how chemical reaction can 
be suppressed for fermionic and bosonic KRb molecules 
under different confinements and electric fields. In sec- 
tion IV, we discuss the importance of inelastic collisions 
of molecules in different confining states. Finally, we con- 
clude in section V. 

In the following, quantities are expressed in S.I. units, 
unless explicitly stated otherwise. Atomic units (a.u.) 
are obtained by setting h = Atteq = 1. 
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II. THEORETICAL FORMALISM 

In this section, we explain the theoretical formalism 
we use. Former studies have dealt with collisions in two 

dimensions lima mm 

but were restricted to small 
confinements or assumed no transitions between confin- 
ing states. In the present formalism, we have no such 
restrictions. Our method is based on a frame transfor- 
mation between spherical to cylindrical coordinates, sim- 
ilar to that employed in Ref. (2|| [26| for example. The 
frame transformation has the advantage of treating in full 
detail the microscopic physics of the molecule-molecule 
interaction, while projecting onto appropriate two dimen- 
sional scattering states. We consider two ultracold polar 
molecules of masses m\,m2 and positions fi,f2 from a 
fixed arbitrary origin O (see Fig.[T]-a). The molecules are 
confined in a harmonic oscillator trap V^ = m T uj 2 z 2 /2 
for molecule r = 1,2, of angular frequency u = 2-kv. 
An electric field applied along the confinement direction 
z polarizes the molecules, giving them dipole moments 
d T = d T z. We use cartesian coordinates [x Tl y T , z T ) to 
describe the vector r T . We also use the center-of-mass 
(CM) coordinate R = (miri + 7712^2 )/{m\ + 777,2) and the 
relative coordinate r = r% — r\ (see Fig. [T]-a). We use the 
cartesian coordinate (X, Y, Z) to describe the vector R, 
and either cylindrical coordinates (p, z, p) or spherical co- 
ordinates (r, 9, ip) to describe the vector r (see Fig. 03-b), 
with p = r sin 9 and z — r cos 9. Both the electric field 
and the harmonic oscillator potential are applied along 
the z axis, which we take as the quantization axis. 

A. Hamiltonian 

The total Hamiltonian of the system is 

fltat = Ti + Tz + V (1) 

with T T = —h 2 'S/ 2 /(2m r ) representing the kinetic energy 
operator of the molecule r. V, the potential energy, is 
given by 

v = v; bs + Kdw + v dd + v; Q =1 + y h ;= 2 

iAe -(r-r mia )/r c _ ^6 d 2 (1 - 3 COS 2 9) 

r 6 47T£o r 3 

+ - (mi uj 2 z\ + m- 2 lo 2 z^j . (2) 

The first term on the right hand side represents an appro- 
priate imaginary potential capturing the overall chemical 
couplings at short-range. It replaces ab initio calcula- 
tions of the electronic structure of trimer and tetramer 
alkali complexes, which remain incomplete for KRb 0- 
[Tol . [27j . For the time being, an absorbing potential 
has shown ver y g ood agreement with experimental re- 
sults [l|, [U, [U |32| for KRb molecules. We use the same 
absorbing potential here. The second term represents 
an isotropic van der Waals interaction ; the third term 
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FIG. 1: (Color online) a) Position vectors of the molecules. 
The electric field is along the z direction, b) Spherical coor- 
dinates (r, 9) and cylindrical coordinates (p, z) of the relative 
coordinate. We suppose tp = in the picture. 



represents the dipole-dipolc interaction where d T repre- 
sents the expectation value in the z direction of the dipole 
moment induced by the electric field ; and the last two 
terms represent the one dimensional harmonic oscillator 
trap that confines the molecules in a plane perpendicular 
to the z direction. The initial energy of the molecule r 
in the trap is given by e„ T = hw(n r + 1/2), where n T 
represents the associated quantum number of the har- 
monic oscillator state they are loaded into. The asso- 
ciated function is the usual normalized eigenfunction of 
the harmonic oscillator g nT (z T ). 



B. Symmetrized internal and external states 

We consider here identical molecules with same masses 
(mi = m 2 ) and same dipoles (d\ = d 2 — d). As the 
molecules are identical, we have to construct an overall 
wavefunction ^ of the system for which the molecular 
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permutation operator P gives 



Appendix A) 



(3) 



with ep = +1 for bosonic molecules and ep = —1 for 
fermionic molecules. This overall wavef unction <]/ is con- 
structed from an internal wavefunction \a\ a 2 ) repre- 
senting the electronic, vibrational, rotational and spin 
degrees of freedom of molecule 1 and 2 respcctivclly 
; from an external wavefunction \n1n2) representing 
the one dimensional individual confining wavefunction 
ffni (^l) <7n 2 (^2) ; and finally from a two dimensional colli- 
sion wavefunction in the plane perpendicular to the con- 
finement. 

We first build symmetrized states of the internal wave- 
function 



\ai a 2 ,rj) 



1 



V2U 



\ai a 2 ) + rj\a 2 oti) 



(4) 



for which P \a\ a 2 ,T}) = f] \ a\ a 2 , r]}. r\ is a good quan- 
tum number and is conserved during the collision. If the 
molecules are in the same molecular internal state, only 
the symmetry 77 = +1 has to be considered. If they are in 
different internal state, both symmetries r\ = ±1 have to 
be considered. We omit explicit reference to the internal 
wavefunctions \aia 2 ,rj) in the following, but the quan- 
tum number 77 still plays a role in the selection rules, as 
discussed in Appendix C. 

We next build symmetrized states of the external con- 
fining wavefunction 



I m ri2,7) = 



f 



^2(1 + 5 nun2 ) 



|ni n 2 ) + "f\n 2 m) 



(5) 



with P\n\ 77,2,7) = 7 \n\ 772,7). 7 is a good quantum 
number and is conserved during the collision. If the 
molecules are in the same external confining state, only 
the symmetry 7 = +1 has to be considered. If they are 
in different external state, both symmetries 7 = ±1 have 
to be considered. It is useful at this point to turn into a 
rclative/CM representation of the confining states. It is 
easy to show that the Hamiltonian (JTJ) can also be written 
in the relative/CM representation as 



T 



T, 



CM 



abs 



vdW 



ho 



CM 
ho 



(0) 



with T re i = -h 2 V 2 P /(2p) and T CM = 
p = mim 2 /(mi + m 2 ) and m to t = 



= -ft 2 V|/(2m tot ) 
mi + m 2 , K r „ o1 ~ 



"ho — 

puj' 2 z 2 /2 and V^ M = m tot uj 2 Z 2 /2. The associ- 
ated energies and functions will be denoted £ n ,£jv and 
9n( z ), 9n(Z). These harmonic oscillator states in the 
relative and CM coordinates are related to those in in- 
dependent particle coordinates g ni (zi), g n2 (z 2 ) by (see 



ffm(zi) 9n 2 (z 2 ) 



v / 2 2 ("i+"2) ni\n 2 \ 

m 712 min(fcjfc ) min(ni — fc,U2 — fc') 

EE E E 



fc=0fc'=0 q=0 



q'=0 



ni\n 2 \ 



{k - q)\ (V - q)\ q\ {n\ — k — q')l q'\ (n 2 



q')\ 



x (_ 1 )«i-fe 2 9 2 ? y/¥n\y/2 N N\g n {z)g N {Z) (7) 



with 



n = —2q' + n\ + n 2 
N = -2q + k + k'. 



k-k' 



(8) 



We give in Appendix A explicit relations between \ni n 2 ) 
and |n, N) states for low values of quantum numbers, < 
n ii n 2 < 2, and in Appendix B the relations between the 
symmetrized individual representation \rii 712,7) and the 
rclative/CM representation \n, N) states, using Eq. ([7]) 
and Eq. ©. 



C. Diabatic-by-sector method 

To solve the Schrodingcr equation for we work in the 
relative/CM representation |n, N) since we know how to 
come back to the physical \ni 712,7) representation. In 
the relative/CM representation, the collisional problem 
depends only on the coordinate Z and the relative vector 
r, and not on the coordinates X and Y. In the follow- 
ing, we explicitly remove these two coordinates from the 
problem. If we use the coordinate Z and spherical coor- 
dinates to represent r, the Hamiltonian is given by 



H 2 1 d 



d 



L' 1 



2p r 2 dr \ dr 



2pr 2 

h 2 



Kibs 



vdW 



ho 



2m tot dZ 2 



CM 



ho 



v dd 



(9) 



If we use the coordinate Z and cylindrical coordinates to 
represent r, the Hamiltonian is given by 



H 



2ji \df + 
2p dz 2 



1 d_ 

P dp 



1 d 2 

p 2 dip 2 

h 2 



+ Kibs + l-vd W + Vdd 



if 



2m to t dZ 2 



+ (io) 



In a diabatic-by-sector method [28l - l30fl . using a spher- 
ical coordinate representation of the wavefunction, the 
range over the Schrodinger equation to be solved, r m i n < 
r < r 
(r, 



max , is divided into N„ sectors of width Ar = 
- r m i R )/N s . The middle of each sector corresponds 



to a grid point r p , with p = 1,...,N S . At each grid 
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point r = r p , we use Ni normalized Legendrc polyno- 
mials P^f L (cos 6>) for a given value of Ml, the quan- 
tum number associated with the azimuthal projection 
of the orbital angular momentum L on the z direction, 
to diagonalize the angular Hamiltonian T-L ML,n (r,9) = 

L 2 /(2pr 2 ) + Vabs + Kdw + Vdd + Vho 1 of thc Hamiltonian 
in Eq. ^ . The resulting eigenfunctions are the adiabatic 
functions xf L ' r> \r p ;9) with j = l,...,Ni. They are used 
as a basis set for the representation of the total wave- 
function 

~ E x]!y\r p -e)g N (Z)F^ N (r p] r) t— (H) 

for a given adiabatic state j. The associated eigenener- 
gies of the angular hamiltonian are the adiabatic energies 
Ej(r p ). They converge to the relative harmonic oscillator 
energies e„ with n = 0,...,Ni — 1 at large r p , so that a 
one-to-one correspondence can be identified between the 
adiabatic quantum states j = l,...,Ni and the relative 
harmonic oscillator quantum states n = 0, TV; — 1. In 
practice, we use a truncated number of adiabatic func- 
tions iVadiab -C iVj . If we restrict the independent oscilla- 
tor quantum numbers < m, n-x < ^osc X i then the maxi- 
mum value that the relative quantum number n can take 
is 2n™f and we choose iV adiab = 2n™. In Eq. ([IT]), 
we use the fact that there are no terms in © that cre- 
ate mixings between different values of N. Moreover, the 
potential V does not depend on the azimuthal angle (p. 
As a consequence, the quantum numbers N and Ml are 
conserved during thc collision. 

The total energy E is equal to e ni + £n 2 + E c , where 
£ni , £n 2 are the energies of thc molecules 1, 2 in the con- 
fining potential, when they start initialy in ni,n%, and E c 
is the initial collision energy between the two molecules 
in the two dimensional plane. E is conserved during 
the collision. Solving the time-independent Schrodinger 
equation = E^ provides the following set of close- 
coupling differential equations in spherical coordinates 
for each values of Ml, t) and N, from a state j to a state 

f 

+ £ Uffp\r p ;r)F^ N (r p ;r) = (12) 
j"=i 

where 
Up,«{r p ;r) = 

r X-' L '\r p ;0)H M ^(r,e) X A J-'\r p ;9) sin6d6. (13) 
Jo 

The goal is to find all the elements Fjtj" v ' ■ We em- 
ploy thc standard method of thc propagation of the log- 



derivative matrix [33| 

Z M LtV ,N {rp . r) = 

Ud/dr)F ML >'i> N (r p \ r)^F AlL ^ N (r p ; r) J (14) 

with matrix elements Zjfj" v,N (r p ;r), and obtain these 
elements for all possible states j to all possible states j'. 
In the diabatic-by-sector method, one has to perform a 
transformation operation from sectors to sectors, since 
the adiabatic functions x ML ' ri { r p\&) change from r p to 
r p+ i. Then the log- derivative expressed in the basis of 
the sector p+ 1 at the distance r = r p + Ar/2 separating 
the sector p and p + 1 is given by 

Z M ^ N (r p+1 ;r = r p + Ar/2) = 

P Z ML ' r i' N (r p ;r = r p + Ar/2) P' 1 (15) 

with the passage matrix 

Pj'j = Xr L ' V (r p+1 ;9) X f L ' n {r p ; 9) sm9d9. (16) 
Jo 

D. Asymptotic matching 

Compared to free molecules in 3D, the external con- 
finement V^® 1 in Eq. ([9]) persists at large intermolecular 
separation r, and the spherical representation of r is not 
appropriate anymore. Instead, we use in the asymptotic 
region cylindrical coordinates appropriate to the poten- 
tial ■ For a given state of relative quantum number 
n, we now expand the total wavefunction as follows 

*^' N (p,z,<p,Z) = 

Y.9n»{z)9N{Z)G^ N {p) e —. (17) 

In the following, we will use the short-hand notation 
£ = Ml,t], 7, N. Note that because we use the coor- 
dinate Z and the wavefunction g^{Z) in both spheri- 
cal and cylindrical representation, the external confine- 
ment V^ M is always well described. When p — > oo, 
Vkbs + Vvdw + Vdd — > 0, and the close-coupling asymptotic 
Schrodinger equations become 

f _ ti^d^ h 2 {Ml -1/4) 
\ ~ 2p,dp^ + 2p~f 

+ e n + e N -E^G{, n {p) = 0. (18) 

At large p, the radial function G^ n , n {p) in Eq. (jTTJ) is a 

linear combination of two possible solutions G n , (p) of 
Eq. (JTHJ) , and takes the form 

Gin (P) -* (p) 5 n , n i + Gf;, (2) (p) Kl n . (19) 
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K^, n represents an element of the reactance matrix. 

The functions G^, represent the regular and irregular 
asympotic solutions of the radial Schrodingcr equation 
Eq. JTT 



G^ip) = p 1/2 JM L (kn>,N P ) 

p 1/2 N Ml {K',n p) 



G^\p) 



(20) 



,N 



where Jm l , Nm l are Bessel functions [38[ and k„ 
\Jl p[E — e n > — £n)/}i is the wave number in the channel 
n' of the relative harmonic oscillator. If E — e n > — £n < 0, 
the modified Bessel functions have to be used instead. 

To determine K, we must transform between the 
spherical wavefunction that captures the short-range 
physics and the cylindrical wavefunction that captures 
the asymptotic boundary conditions. The regular and ir- 
regular spherical radial functions F^'^' 2 \r p ; r) and their 
derivatives can be connected to their cylindrical asymp- 
totic counterpart G^'^' 2 '(p) by equating the wavefunc- 
tions Eq. (TlTj) and Eq. (JTTJ) and their derivatives at a 
constant sphere of radius r = r max 



4f 2) (r>^;r) 



Xj- {r P =N s 



^-g n ( z )G^ 2 \p) S m6d6 



(21) 



l( Ff,f> 2 \r p = Ns ;r) 



dr \ 



j .i 



(22) 



with the one-to-one correspondence {n = 0, iV a diab — 
1} = {j = 1, Afadiab} between the quantum numbers n 
and j. r p= N s is the middle of the last sector N s . This is 
a similar matching procedure that connects short-range 
democratic hyperspherical coordinates to asymptotic Ja- 
cobi coordinates employed in atom-molecule chemical re- 
active scattering studies (2^-[30j. Convergence with re- 



spect to N; 



adiab 



and 



is found when the Wronskian 



matrix 



dr \ 



or 



(23) 



converges to the unit matrix. 

The K matrix is determined by the matrix operation 



Z*Ft<W - (d/dr)(F^) 



(24) 



The scattering matrix S in the relative/CM representa- 
tion is determined by 



(25) 



where in this equation, I represents the unit matrix. The 
scattering matrix in the symmetrized individual repre- 
sentation |niTi2,7) is found by gathering all individual 
scattering matrices S corresponding to different values 
of N and by applying a transformation from the rela- 
tive/CM representation to the symmetrized individual 
representation 



gM L ,vn = U i S Ml ' 71 '~ < ' N X 



U 1 



(26) 



The transformation matrix U, with elements 
,N = (nin 2 ,^\n,N), can be found using 
the relations in Appendix B. We use the transpose U T 
of the matrix U instead of its inverse because U is not 
generally a square matrix. 



E. Observables 

After a collision, the quantum probability from an ini- 
tial state ni n-2 to a final state n' 2 for defined numbers 
M L , 77 , 7 is given by P* L >™ = \S M , L '? n I 2 . The 
clastic, inelastic (confining state changing) and reactive 
probabilities are given by 



pel, Ml ,77,7 _ pMhtVtl 

ni n 2 ,ni n 2 



pin, Ml,t\,1 
pic, Mi,,?;, 7 



Y _ p ci , Ml ,V,1 



pML,V,1 



n\ n ,n\ n 2 

pin, M L , n,7 



(27) 



We mean by "inelastic", processes that change the ex- 
ternal confining states of the molecules. Finally, for an 
initial state n-y n 2 , the elastic, inelastic, and reactive cross 
sections are given by j33 - [3r| 

V 11 - qMLV'-r I 2 x A 

/ / I 1 D n 1 n 2 ,n 1 n 2 \ X ^ 
Ml,V,1 

pin, M L , 11,1 x ^ 

Ml,ti,1 

p.*w.7 X A. (28) 

Ml,V,~I 

The inelastic state-to-state cross section is given by 

h 

W, 

Ml,V,1 



h 
h 

V^pEc 



n± ri2 to n' 2 



_L = y p M L ,n,-i xA . 



(29) 



The factor A represents symmetrization requirements for 
indistinguishable particles in a same internal and con- 
fining state [III H3- The cross sections are found by 
summing over all the contributions of different values 
of M^, 77, 7. For the ultralow energies involved in this 
study, only the first partial wave will be required for in- 
distinguishable molecules (same internal states r\ = +1 
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and same confining state 7 = +1): the Ml = par- 
tial wave for indistinguishable bosons and the Ml = ±1 
partial wave for indistinguishable fcrmions. The temper- 
ature dependence of the loss rates in the two dimensional 
plane is found by averaging the cross sections over a two- 
dimensional Maxwc ll-Boltz mann distribution of the rel- 
ative velocity v = y/2E c /[j, in the two-dimensional plane. 
This gives a two dimensional thermalized rate 

/>oo 

^tf' re = / <tni C vf(v)dv (30) 
JQ 

with 

/(i^-^-oe" 5 ^ (31) 
3 y ' k B T y ' 

where kp is the Boltzmann constant. The rate in Eq. 1301 
corresponds to the rate per molecule, not the event or 
collision rate [i~lll37j. 

Selection rules apply due to symmetrization of the 
wavefunction under permutation of identical molecules 
(Appendix C). The rules are 

V1(-1) Ml = V (-1) L = V (-l) ML+n = e P . (32) 

This limits the summation over Ml, 77, 7 in Eq. (|28|) 
and (f30|) and the values of the quantum numbers j" and 
n" used in Eq. ((TTJ) and Eq. (fT7|) . 

In the following, we will consider molecules of KRb as 
an illustrative example of experimental interest [3, 0, 
[22| . For concreteness, we will take the isotope 39 K 87 Rb 
for the bosonic molecules ; the results for the bosonic iso- 
tope 41 K 87 Rb @ are nearly identical. We take the iso- 
tope 40 K 87 Rb for the fermionic molecules (HHE^. Con- 
vergence of the results have been checked with the match- 
ing distance r max and the number of adiabatic functions 
-Wadiab included in the expansion of the wavefunction. 
Unless stated otherwise, we choose r m ; n = 10 ao and 
r max = 10000 ao (ao — 0.529 Angstroms is the Bohr 
radius), N s = 10000 sectors, < ni,n 2 < rCf = 3, 
^Vadiab = 2 n™ x = 6 and only the first partial waves 
Ml = 0,1 depending on the species and the selection 
rules involved. We used Ni = 80 Lcgcndrc polynomials 
for v < 100 kHz and N t = 120 for v > 100 kHz, to 
construct the adiabatic functions. This yields converged 
results of 10 % at most for the elastic rates (more es- 
pecially at high confinement) and 1 % for the reactive 
and inelastic rates. For Vabs, we use A = —10 K and 
r c = 10 ao, which adequately reproduces experimental 
loss rates in three dimensional collisions [2f . 



III. SUPPRESSION OF CHEMICAL 
REACTIONS 

We discuss in this section how chemical reactions pro- 
ceed when the reactants are subject to different confine- 
ments and electric fields. We present in Fig. the adi- 
abatic energies £j(r p ) for the symmetry j (— 1) Ml = — 1 



(upper panel) and the symmetry 7 (— 1) Ml = +1 (lower 
panel), for a trap with v = 20 kHz and induced dipole 
moment d = 0.1 D. These energies converge at large r 
to the energies of the relative harmonic oscillator e n . To 
associate a specific confined collision with a symmetry 
j(—1) Ml , one has to use Eq. (f3"2")) . If the molecules are 
identical fcrmions in the same internal state, r\ — +1 and 
€p = —1, and then 7 (— 1) Ml = —1, so the scattering 
problem only employs the black and red dashed curves 
of the upper panel in Fig. O In addition, if the identical 
fermionic molecules are in the same external state, then 
7 = +1, and the scattering problem only uses the black 
curves. If however the identical fermionic molecules are 
in different internal states, both values of 77 are relevant. 
Then, in the case of 77 = — 1, now 7 (— 1) AIl = +1, and 
the black and red dashed curves of the lower panel have 
to be employed as well. If the fermionic molecules are in 
different internal states, but in the same external state, 
then 7 = +1, and one has to use only the black curves of 
both panels. 

Using similar arguments, if molecules are identical 
bosons in the same internal state, one has to use the 
black and red dashed curves of the lower panel. If 
besides they are in the same external state, only the 
black curves have to be used. If they are in different 
internal states, all black and red dashed curves of both 
panels have to be used, while only the black curves 
of both panels are used if the identical bosons are in 
different internal states but in the same external state. 
The case of two different polar molecules corresponds to 
all curves of all symmetries employed. Also, note that 
because 7(-1) a/l = (-1) L = (-l) M ^+ n in Eq. (j32|l. 
the values of L and Ml + n are odd for the upper 
panel and even for the lower panel, and the 7 = +1 
(7 = —1) curves corresponds to even (odd) relative 
quantum numbers n (7 = (—1)"). Therefore, symmetry 
consideration are essential for the dynamics of ultracold 
molecules in confined geometry and electric field. 

We now discuss the differences between the symme- 
tries rather than a specific confined collisional case. We 
focus on the symmetry 7 (— 1) Ml = —1 with 7 = +1 
(black curves of the upper panel in Fig. [2]) and on the 
symmetry 7 (— 1) Ml = +1 with 7 = +1 (black curves 
of the lower panel in Fig. [2]). The former case corre- 
sponds to the dynamics of identical indistinguishable 
fcrmions and the latter to the dynamics of identical 
indistinguishable bosons. By indistinguishable, we mean 
identical molecules in the same internal and external 
states. For the discussion, we focus only on the lowest 
black curve if we assume molecules in the ground state 
of the trapping potential. Two striking differences can 
be seen due to the statistics of the systems. First, the 
lowest curve connects at short distance to an adiabatic 
curve with a L = 1 adiabatic barrier V& (depicted 
with a green arrow) for the 7 (— 1) AIl = — 1 symmetry 
(upper panel), while no barrier is present (L — 0) for 
the 7 (— 1) AIl = +1 symmetry (lower panel). This 
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FIG. 2: (Color online) Adiabatic energies versus r for 
the 7 (— 1) Ml — — 1 symmetry (upper panel) and for the 
7( — 1) ^ — +1 symmetry (lower panel), for v — 20 kHz 
and d = 0.1 D. The black (red dashed) curves correspond to 
7 = +1 (7 = —1) manifolds. We also show how values of L 
and n adiabatically connect. 14 is the height of the barrier 
for molecules in the lowest confining state (n = 0). 



indist. fermions 




3-1 ^ 



ndist. bosons 




makes indistinguishable bosonic molecules likely to 
chemicaly react in confined geometry compared to 
fermionic molecules. Second, the lowest curve (7 = +1) 
corresponds to Ml = ±1 for the first symmetry while 
it corresponds to Ml = for the second one. Under an 
electric field, the Ml = component always corresponds 
to an attractive dipole-dipole interaction whereas the 
Ml = 1 component corresponds to a repulsive dipole- 
dipole interaction (which can eventually turns into an 
attractive one at higher dipoles @, [ll| ) . For this rather 
small confinement, it means that we can still, up to a 
certain dipole, use an electric field to increase the barrier 
Vi, for indistinguishable fermions. This is not true for 
indistinguishable bosons. We will refer to this kind of 
suppression as "statistical suppression", as it depends 
on the fermionic/bosonic character. To get suppression 
for indistinguishable bosons, we will have to increase the 
confinement and the electric field, which will be refered 



FIG. 3: (Color online) Height of the adiabatic barrier V& ver- 
sus d and v for indistinguishable fermions (upper panel) and 
for indistinguishable bosons (lower panel) in the lowest con- 
fining state. 



in the following to as "potential energy suppression" . 

To understand these two types of suppression, it 
is useful to plot the height of the barrier Vb, which 
the molecules at ultralow temperature must tunnel 
through. We plot this barrier in Fig. [3] for the symmetry 
7(— 1) Ml = —1 with 7 = +1 (upper panel) and for the 
symmetry 7 (— 1) Ml = +1 with 7 = +1 (lower panel), as 
a function of the confinement v and the dipole moment 
d induced by the electric field, for the lowest confining 
state. For the first symmetry (upper panel), there are 
two ways to get a high barrier. One way is for small 
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indist. fermions 




indist. bosons 




v m4 10 

FIG. 4: (Color online) Elastic and reactive rate coefficient 
versus d and v for indistinguishable fermions (upper panel) 
and for for indistinguishable bosons (lower panel) at E c = 
500 nK. The elastic curve is plotted in red. 



confinements and small d. The barrier increases to reach 
a maximum at d ~ 0.15 D. The fact that the barrier 
decreases for higher dipoles comes from contributions of 
higher values of L = 3, 5, ... @, [10]. For d w 0.15 D, if 
we follow this maximum of the three-dimensional plot for 
increasing confinements, we see that Vb decreases again. 
When v increases, the zero-point energies (the ones at 
large r in Fig. [5]) increase while the barrier is not affected 
at short distance because the confinement is small. Then, 
the effective height of the barrier is decreased [l5| as v 
increases. The second way to achieve high barriers Vj, 
is for high dipoles and high confinements. The barrier 



increases monotically, emphasizing the electric field sup- 
pression of confined chemical rates. When the molecules 
are highly confined in a two dimensional plane perpendic- 
ular to an applied electric field, they collide side-by-sidc. 
This repulsive electric interaction enhances the barrier 
and makes the molecules stable against collisions (l2l - fl7| . 

For the second symmetry (lower panel), there is only 
one way to increase the barrier. The striking difference 
is that for small confinement and/or small dipoles, there 
is no barrier at short range as already seen in Fig. [2] The 
only way to raise the barrier is for high confinements and 
high dipoles as for the first symmetry, where the electric 
dipole repulsion come into play. The rise of the barrier 
at high confinements and high dipoles is independent of 
the symmetrization of the molecules, as Vb converges to 
similar values for both cases. 

The behavior of Vb has crucial consequences on the dy- 
namics of the molecules. To get the rate coefficients of 
a specific confined collision, one has to add the rates ob- 
tained from a scattering calculation using the adiabatic 
curves of the individual symmetries 7 (— 1) Ml involved 
in the specific problem. The rates for the symmetry 
7 (— 1) Ml = — 1 with 7 = +1 is presented in the up- 
per panel and for the symmetry 7 (— = +1 with 
7 = +1 in the lower panel of Fig. [4] as a function of v 
and d for a collision energy E c = 500 nK. Qualitatively, 
the behavior of the reactive rates is opposite to the height 
of the corresponding barriers, while the elastic rates in- 
crease only in a monotonic way with d and v. For small 
confinements and dipoles (say v = 20 kHz, d = 0.15 D), 
the reactive rates are suppressed for the first symmetry 
(upper panel) representing approximatelly 10~ 2 of the 
elastic rates. No such suppression is seen for the second 
symmetry (lower panel). This shows that this statistical 
suppression is only due to symmetrization requirements, 
but has the advantage to work at rather realistic experi- 
mental confinements and dipoles. For high confinements 
and dipoles, the reactive rates of fermions and bosons can 
be suppressed by three to four orders of magnitude com- 
pared to the ones at small confinements. This is made 
possible by the anisotropy of the dipolar interaction of 
polar molecules in confined geometries as explained in 
Refs. [H-[l7|. 

The elastic rates increase as d 4 or d, depending on the 
collision energy and magnitude of the dipole [39(, and 
increase with v [l5|, Therefore, this potential en- 

ergy suppression of the reactive rates and enhancement 
of the elastic processes will help evaporative cooling of 
fermionic and bosonic molecules, and will make amenable 
the creation of degenerate Fermi gases or Bose-Einstein 
condensates of polar molecules. This suppression is not 
due to symmetrization requirements but to the fact that 
the molecules possess a permanent electric dipole mo- 
ment. Therefore, this suppression will also be effective 
for molecules in distinguishable states or even for non- 
identical polar molecules. 

It is worth noting that the fermionic statistical sup- 
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FIG. 5: (Color online) Rate coefficient versus collision 

energy E c for d = 0.1 D and v — 20 kHz, for indistinguishable 
fermions (upper panel) and indistinguishable bosons (lower 
panel) , initially in the ground state of the confining trap n\ = 
ri2 = 0. The thick solid (dashed) curve corresponds to reactive 
(inelastic) scattering. The thin solid black lines represent the 
confining state-to-state rate coefficients. 
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prcssion is still effective if the fermions are in different 
external states (7 = ±1), since both black and red dashed 
curves of the upper panel in Fig. [5] have to be used. The 
red curves corresponds to a Ml = component, whose 
barrier height VJ, decreases for increasing electric field. 
There is no statistical suppression at all if the molecules 
are in different internal states (77 = ±1), because the 
curves from the lower panel in Fig. [2] have to be used 
including the barierless curve L = 0. This has been con- 
firmed experimentally (22| . 

Finally, no statistical suppression can occur in the case 
of different polar molecules, for which all curves of all 
symmetries in Fig. [2] should be employed. Only the po- 
tential energy suppression can apply in that case. 



FIG. 6: (Color online) Thermalized rate coefficient versus 
d for T = 500 nK and v = 20 kHz. The solid (dashed) 
curves correspond to reactive (inelastic) processes. The red 
(blue) curves correspond to fermions (bosons) in same in- 
ternal states, but not necessarily in same external states. 
The molecules are considered initially in n\ = 0, ni — 
(upper panel), in ni = l,ri2 = 1 (middle panel), and in 
n% = 0,n2 = 1 (lower panel). 



IV. INELASTIC COLLISIONS BETWEEN 
CONFINING STATES 

We saw that chemical suppression of indistinguishable 
fermions and bosons can always be obtained if sum- 
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cicntly high confinements and electric fields are applied. 
However, the magnitude of these high confinements is 
still beyond of those that can be currently achieved 
experimentally. For a realistic experimental frequency 
of v ~ 20 kHz, loss of indistinguishable fermions can 
be suppressed taking advantage of the alternative 
statistical suppression whereas loss of indistinguishable 
bosons cannot realistically be suppressed. Moreover, 
for small confinements, it is possible that higher trap 
confining states can be populated. The reason is that 
the energy spacing between two allowed confining states 
As = 0.96 /J.K for v ~ 20 kHz can be of the order of 
the temperature T ~ 500 nK of the gas. Then, the 
changing-state dynamics of molecules in small confining 
optical lattices must be understood as well. We consider 
in the following fermions and bosons in same internal 
states but not necessarily in the same external confining 
states, for a realistic confinement of v = 20 kHz. 

We present in Fig. [5] the non-thermalized rate coeffi- 
cients /3oo' le = °oo' re v as a function of the collision en- 
ergy for the inelastic and reactive processes, for indistin- 
guishable fcrmionic molecules (upper panel) and indis- 
tinguishable bosonic molecules (lower panel) in same in- 
ternal and external states. The molecules start in n,\ = 
and ri2 = and v = 20 kHz, d = 0.1 D. States be- 
tween < 711,712 < n™ x = ^ have been used for 
collision energy E c > 1 yuK to converge these results. 
At ultralow energy, the fermionic reactive rate scales as 
E c , and as ln~ 2 (-y/2 \i E c ) for the bosons, in agreement 
with the threshold laws |4l[. When the collision 
energy is sufficiently high, excited confining states be- 
come energetically open. Overall, bosons react at higher 
rate than fermions, as expected, since there is no bar- 
rier for bosons, whereas there is a barrier for fermions. 
Moreover, molecules that start in the ground confining 
state are much more likely to react chemically than to 
go to higher confining state. The inelastic rate for the 
fermionic molecules is an order of magnitude smaller than 
its reactive rate. It is a factor of 3 — 8 smaller than the 
reactive rate for the bosonic molecules. 



A. Gas in thermal equilibrium 

We now consider a thermal equilibrium at a tempera- 
ture of T = 500 nK. The population p of the molecules 
in n T is given by a Maxwell-Boltzmann distribution 

e k B T 

p(n T ) = ZIEZ- ( 33 ) 

V e fc s T 

At T = 500 nK in a trap with v = 20 kHz, 
p(n T = 0) ~ 0.852, p(n T = 1) ~ 0.126 and 
p(n T = 2) ~ 0.019. In the following we will ne- 
glect contribution of molecules in n T = 2, and consider 
only molecules in n T = 0,1 for simplification. The 
coefficients p(n T ) will play a role in the rate equations 



below. 

We present in Fig. [6] the thermalized rates Pqq C and 
Pqq 111 (upper panel), P^{ re and f3i{ ln (middle panel), and 
/3j{ re (lower panel), as a function of d for v = 20 kHz at 
T = 500 nK. The reactive and inelastic rates are plot- 
ted as a thick and dashed solid line. The fermionic and 
bosonic case are plotted in red and blue respectively. 

We discuss first the case of molecules in the ground 
states ni = 0, ri2 = (upper panel). For bosons, the reac- 
tive rate is high and the inelastic collision is insignificant. 
For fermions however, the inelastic rate can reach 20 % of 
the amount of the reactive rate at d = 0.23 D. The mag- 
nitude of the thermalized inelastic rates is proportional 
to the amount of molecules allowed by the Maxwell- 
Boltzmann distribution at T = 500 nK to have kinetic 
energy greater than the first excited inelastic thresholds 
ni = 1.^2 = 1 and n\ = 0, no, = 2 at 1.92 /iK. We 
also plot in circles (fermions) and triangles (bosons) the 
non-thermalized reactive rate /3qq = (Tqq v. Wc sec that 
Pqq 8 = /3qq is a reasonable approximation at small dipolc 
moments. /3qq differs by 35 % from /3^J 1C at the high- 
est dipole, however. This comes from the fact that at 
these dipoles, the molecules do not collide in the Wigner 
regime anymore and the height of the barrier for fermions 
(or characteristic energy for bosons) is comparable to the 
temperature. Note that if the confinement is increased 
to v = 30 kHz, the inelastic rate (represented as thin 
dashed black lines) decreases by about an order of mag- 
nitude, because for a same temperature, it is harder to 
excite molecules in higher confining states as the energy 
thresholds increases with the confinement. Then the in- 
elastic collisions for ground state molecules become less 
important as the confinement increases. 

If the molecules are now in the first excited states 
n\ = 1,H2 = 1 (middle panel), reactive collisions, for 
both bosons and fermions, are about 30 % smaller than 
the ones for molecules in n\ = 0, n% = 0. A qualita- 
tive explanation is that ni = l,n2 = 1 (which has a 
7 = +1 symmetry) projects onto an n — 0, ./V = 2 state 
and a n = 2,N — state (see Appendix B). When we 
look at the corresponding adiabatic energies in Fig.[2]for 
the fermions, the n = 2 curve connects to the L = 3 
adiabatic barrier which is much higher than the L = 1 
barrier, suppressing more strongly the reactive collisions 
and increasing inelastic collisions. For bosons, the reac- 
tive rates are still high compared to fermions, because 
the n = curve connects to a L = curve. However, 
the reactive rates are smaller than for the n\ = 0, ri2 = 
case, because there is now the n = 2 curve that connects 
to a L — 2 curve, suppressing chemical reactivity. The in- 
elastic processes are much important in the present case 
because the Maxwell-Boltzmann distribution allows all 
molecules to have sufficient kinetic energy to contribute 
to the inelastic process, while in the precedent case, only 
a part of the molecules were allowed to contribute to the 
inelastic process. For bosons, the inelastic magnitude is 
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about half the reactive rate (at most, at d = 0.3 D), 
but for fermions, it can even exceed the reactive rate for 
d > 0.1 D. 

Finally, we discuss the case of molecules in different 
states n\ = 0, ri2 = 1 (lower panel). This channel cannot 
decay to the energetically allowed n\ = 0, n-2 = channel, 
because the two channels correspond to different values 
of N. However, the molecules arc in different confining 
states now so that two contributions 7 = ±1 are involved 
in the calculation, and both black and red dashed curves 
of Fig.[5]have to be used. This is shown in the lower panel 
of Fig. [6] as thin solid line for 7 = +1 and thin dashed 
line for 7 = — 1. Compared to fermionic molecules in 
the same confining states, the reactive rates are bigger. 
This comes mainly from the 7 = — 1 contribution, which 
corresponds to Ml = head-to-tail attractive dipolar 
interactions (see Tab. UJ. For bosonic molecules in dif- 
ferent confining states, the reactive rates are similar to 
those for molecules in same confining states, except that 
the 7 = — 1 contribution gives an enhancement at high 
dipoles due to the Ml = 1 component of the L = 2 adi- 
abatic curve (see Tab. Hj). This component corresponds 
to an attractive dipolar interaction (see Eq. 8 and Eq. 9 
of Ref. [ll|) and can enhance the reactive rate at high 
dipoles. The L = 2 barrier is high at small dipoles (see 
Fig. [5]) and suppresses the reactive rates. However, the 
strong dependence of d 4< - L+1 ^ 2 ^ of the rates [Tl| leads to 
a d 10 dependence, as shown in the figure, and eventually 
makes a significant contribution at high dipoles. 

We saw on one hand that inelastic processes can be 
important for molecules initially in excited confining 
states, especially for fermions, and that on the other hand 
molecules can chemically react at high rates for molecules 
initially in different confining states, even for fermions be- 
cause they are not indistinguishable anymore. What are 
the consequences of this for the dynamics of a molecular 
gas? This is what we answer in the next subsection. 



B. Rate equations 

The rate equations for the density of molecules n„ x (t) 
in state n T as a function of time are given by 

h (t) = -Po£ e nl(t) -/3^ re n (t)ni(t) 

- 0Mtoii°o(*) + j8utooo*i(t) 
m (t) = -plr n? (f ) - $F no (t) m (t) 

" 0u to 00 to 11 *»<>(*) (34) 

where no(i) (ni(£)) arc the individual densities of 
molecules in state n T = (n T = 1). Similar equa- 
tions hold for n T > 2, but for simplicity, to avoid ad- 
ditional inelastic terms in the equations, we assumed 

Pn-r>2 < Pn T =0,l- 

If we assume a gas in thermal equilibrium for each 
time t, the Maxwcll-Boltzmann distribution implies that 
n (i) = p(0)n to t(i) and ni(i) = n tot (i) (we assume 



p(0) + p(l) ~ 1 in our example), where n to t(*) is the 
density of the total molecules. Then by summing the 
equations above, we obtain the rate equation for ntot(i) 

htot(i) = -|/(0)^ re +/(l)^ii re 

+ 2p(0)p(l)/3 T i ro |n t 2 ot W. (35) 

Inelastic rates cancel each other in the full equation, be- 
cause two molecules go back and forth in m = 0, n 2 = 
and ni = 1,71a = 1, without participating in the loss 
process. Although inelastic collisions are responsible for 
the evolution of the individual density of molecules n (t) 
and ni(i) , they are not responsible for the evolution of 
the total density of molecules in the thermal gas. 

At T = 500 nK, pff* ~ /3^ re but p 2 (l) < p 2 (0) 
so that the second term on the right hand side of the 
equation above can be neglected. As a result the den- 
sity of the total molecules will show a faster decay due 
to a fast rate 2 p(0) p(l) (3q{ 1 ° and a slow decay due to 
a slow rate p 2 (0)/3g^ re . For example for fermionic KRb 
at d = 0.2 D, 2p(0)p(l)/3 r 1 ' rc ~ 1.4HT 6 cm 2 s" 1 and 
P 2 (0)(3qq C ~ 5.10~ 8 cm 2 s _1 . The fast and slow decays 
are due to high inter-states reactive rates (collisions be- 
tween different confining states) and low intra-states re- 
active states (collisions between same confining states). 
The two types of decay can be tuned by changing the 
relative populations p(0) and p(l), by changing the tem- 
perature T and/or the confinement v. Note that even 
if the population of the molecules in different confining 
states are not given by a Maxwell-Boltzmann distribu- 
tion, say for example p(0) = 0.5 and p(l) = 0.5, and is 
independent of time, inelastic rates still cancel each other 
in the equation for the total density of molecules. Again, 
inelastic collisions play a role in the loss of molecules from 
individual trap levels, but do not for the loss of the total 
molecules. These theoretical findings well support recent 
experimental data of confined fermionic KRb molecules 
in electric fields [22| . 

V. CONCLUSION 

We have developed in detail a rigorous time- 
independent quantum formalism to describe the dynam- 
ics of particles with permanent electric dipole moments 
in a confined geometry, by treating the reactive chemistry 
using an absorbing potential. Elastic, reactive and inelas- 
tic rate coefficients can be computed for a given collision 
energy, temperature, confinement and dipole moment 
(or electric field), for a system of fermionic or bosonic 
molecules. The selection rules play an important role for 
the dynamics of confined molecules and have dramatic ef- 
fects on the collisional properties. Different rates are ob- 
tained for fcrmionic/bosonic molecules in same/different 
confining states. Two kinds of suppression can occur for 
chemical reactions: a statistical suppression applies only 
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for fermions at rather small induced dipolcs and con- 
finements realistically accessible in an experiment, and 
a potential energy suppression applies for both fermions 
and bosons at rather high induced dipoles and confine- 
ments. Inelastic rates can be important, even as high 
as reactive rates for molecules initially in excited states. 
However, the inelastic rates do no play a role in the loss 
process of the total number of molecules in a gas, since 
molecules are inelastically excited and relaxed, back and 
forth. Only reactive rates are responsible for the evo- 
lution of the loss of the total molecules. Fast and slow 
decays of the molecules can be seen due to inter-state and 
intra-state confined collisions. This work has been highly 
motivated by recent experiments of KRb molecules in 
confined geometry and electric field, and has proved very 
good theoretical support for the experimental observa- 
tions m. 
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Appendix A: Relation between |m 722} and \n,N) 
In Eq. ((7|), we use the following characteristics [42| 



9n T (x) 



1 



2™- n T \ \ nh 



1/4 



e"^fl"„ T (vWi/Si) (36) 



H nT (x + y) = 2-^/2 £ 



k=0 



k\(n T -k)\ 



H„ (x) H m ^ (x) 



H k (xV2)H nT - k {yV2) (37) 



min (n^ ,m T ) 

E 

fc=0 

^- T T7H_ 2k+mT+nT (x)2 k kl. (38) 



k\{m T - k)\ k\(n T - k)\ 
The individual Ini 712) states are written in terms of 



the relative/ CM \n,N) states by 
[00) 
101) 



10) 
02) 
20) 
11) 



|0,0) 

1 

71 

1 



0,1 



1 



1,0) 
1,0) 



i|0,2) + i=|l,l) + i|2,0) 
i|0,2)--L|l,l> + i|2,0) 
>,2)-^|2,0). 



(39) 



Appendix B: Relation between |m 112,7) and \n,N) 

Using Eq. ^ and Appendix A, the symmetrized indi- 
vidual \ni 71-2,7) states are written in terms of the rela- 
tive/CM \n, N) states by 



00,7 = +l> = |0,0) 

01,7 = +1> = [0, 1) 

02,7 = +l> = -^|0,2) + i=|2,0) 

11,7 = +1) = -^|0,2)-i=|2,0) 

12,7 = +1) = \/3|0,3>-J|2,l) 



|22,7 = 


+1) = 


vfi°. 4 >- 


i|2,2) + y||4,0) 


|01,7 = 


-1) = 


[1,0) 




|02,7 = 


-1) = 


IM) 




|12,7 = 


-1) = 


£|1,2>- V 


1|3,0>. 


Note that (- 




= (— l) n + N 





(40) 



fermions 


V 


L 


7 


Ml 


n 




+ 1 


1,3,5 ... 


+ 1 


1,3,5 ... 


0,2,4 ... 






1,3,5 ... 


-1 


0,2,4 ... 


1,3,5 ... 




-1 


0,2,4 ... 


+ 1 


0,2,4 ... 


0,2,4 ... 






2,4,6 ... 


-1 


1,3,5 ... 


1,3,5 ... 


bosons 


V 


L 


7 


M L 


n 




+ 1 


0,2,4 ... 


+ 1 


0,2,4 ... 


0,2,4 ... 






2,4,6 ... 


-1 


1,3,5 ... 


1,3,5 ... 




-1 


1,3,5 ... 


+ 1 


1,3,5 ... 


0,2,4 ... 






1,3,5 ... 


-1 


0,2,4 ... 


1,3,5 ... 



TABLE I: Selection rules for the dynamics of identical bosons 
and fermions in confined two dimensional geometry. 
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Appendix C: Selection rules 

For initial states n\,ni and final states n'^n^, since 
components of different N do not mix together in the 
collision process, we have 

(_l)™i+™2 = r_i)™!+"2 (4i) 

after a collision. 

At long range, in cylindrical coordinates, if we use 
the symmetrized individual representation \n\ ^2,7), the 
permutation P requires the substitutions Z\ — > z 2l z 2 — > 
Z\ , ip — > if + 7r which leads to the selection rule 

?77(-l) Mi =e P . (42) 



If we use the relative representation |n, N) states, then 
the permutation P requires the substitutions z — > 
— z, ip — > ip + 7r which leads to 

V(-l) ML+n = e P , (43) 

from the properties of the g n {z) functions. At short 
range, in spherical coordinates, using the Legendre poly- 
nomials, the permutation P requires the substitutions 
>7r — 6, ip ^ ip + tv which leads to 

V (-l) L = e P . (44) 

We summarize in Tab. [I] the different selection rules for 
identical bosons and fermions. 
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